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ABSTRACT 

Using volume- limited samples drawn from The Sloan Digital Sky Survey Data Release 
7 (SDSS DR7), we measure the tidal environment of galaxies, which we characterize 
by the ellipticity e of the potential field calculated from the smoothed spatial number 
density 1 -\- S of galaxies. We analyze if galaxy properties, including color, Dn4000, 
concentration and size correlate with e, in addition to depending on 1 -\- S. We find 
that there exists a transition smoothing scale at which correlations/anti-correlations 
with e reverse. This transition scale is well represented by the distance to the 3rd- 
nearest-neighbor of a galaxy in a volume limited sample with Mr < — 20 which has 
a distribution peaked at ~ 2 h~^Mpc. We further demonstrate that this scale corre- 
sponds to that where the correlation between the color of galaxies and environmental 
density 1 + ^ is the strongest. For this optimal smoothing Rq no additional correlations 
with e are observed. The apparent dependence on tidal ellipticity e at other smoothing 
scales Rs can be viewed as a geometric effect, arising from the cross correlation be- 
tween {l-\-So) and e{Rs). We perform the same analysis on numerical simulations with 
semi- analytical modeling (SAM) of galaxy formation. The e dependence of the galaxy 
properties shows similar behavior to that in the SDSS, although the color-density cor- 
relation is significantly stronger in the SAM. The 'optimal adaptive smoothing scale' 
in the SAM is also closely related to the distance to the 3rd-nearest-neighbor of a 
galaxy, and its characteristic value is consistent with, albeit slightly smaller than for 
SDSS. 
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1 INTRODUCTION 

While the standard picture of the hierarchical build-up of 
large-scale structure in the universe has been supported by 
ever increasing observational evidence, it is still an ongoing 
task to understand thoroughly the detailed baryonic physics 
associated with galaxy formation, such as gas accretion, star 
formation and feedback. Direct observation of these gaseous 
processes is challenging. On the other hand, given the large 
galaxy samples currently available to us, such as SDSS DR7 



it is the star- format ion related properties, such as color and 
emission-line flux, that are more directly correlated with en- 
vironmental density (Kauffmann et aL]|2004 Blanton et al 



2005b Christlein & Zabludoff 2005). At fixed color, the 



residual dependence of morphology on the environmental 



density is rather weak (Park et al. 2007 Ball et al. 2008 



mann et al. (2004), Blanton k Berlind (2007) and Park et 



Bamford et al.||2008). It is also emphasized by, e.g., Kaufi'- 



al. (2007), that the environment dependence is in fact quite 



( Abazajian et al.||200 9), statistical analysis of galaxy prop- 
erties has proved to be invaluable to probe the underlying 
physics affecting galaxy formation and evolution. 

The morphology-density correlation for galaxies, un- 
covered already in the 1930s ( Hubble||1936 ), shows a clear 
segregation, with early-type galaxies preferentially existing 
in high density regions (^Oemler^^l974^ ^Dressler^^l978| [Post- 



local. While galaxies with abundant close neighbors have 
clearly different properties from those of isolated galaxies, 
on scales larger than ^ lh~^ Mpc, environmental effects 



man fc Geller |[T974l [Blanton et al. ||2005b| [Kauffmann et 



al. 112004 1' Weinm annet al.||2006| [Blanton &: Berlind 1120071 
Park et al. 2007j). Further investigation has indicated that 



probably have little influence on galaxy properties (Blanton 
|fc Berlind |2"007t . 

Theoretically, it has long been known from extensions 
of Press- Schecter theory that massive halos preferentially re- 
side in dense environments ( Bond et al.[[1991 Mo & White 
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1996). Galaxies therein have properties that are different 
from those in small halos. On the other hand, the excursion 
set theory with sharp-k fllter and with dynamics described 
by the spherical collapse model predicts that the formation 
history of halos of fixed mass should depend only on smaller- 
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scale structure not on their large-scale environment (e.g., 
White 1996|). Incorporating the ellipsoidal collapse model 



into excursion set theory gives halo mass functions and bias 



in better agreement with numerical simulations (Bond & 
Myers||1996] |Sheth k Tormen||20Qlt . In its simple version, 



however, the theory takes into account the ellipsoidal col- 
lapse by considering the average tidal field, and thus also 
predicts no correlations between the formation history of 



halos and their specific large-scale environment ( [Sheth fc 
Tormen|200T| 



Recent numerical simulations show that halo proper- 
ties can be strongly correlated with their environment. Halo 
assembly bias reveals that old halos tend to cluster more 
strongly than their younger counterparts, and this bias is 
particularly significant for low mass halos (|Gao et aL 



2005 



Harker et al . 2006; W echsler et al.|2006|[Gao &: White|2007 
Croton et al.,2007, |Jing et al.|2007| ). Using filters other than 
the sharp-k filter in the excursion set theory leads to non- 
Mar kovian random walks, which in turn generate correla- 
tions between halo formation history and large-scale envi- 
ronment (e.g.. Ma et al.|2010 ). Desjacques (2008) shows that 
under the spherical collapse model, the correlations aris- 
ing purely from this non-Markovianity are expected to be 
stronger for more massive halos, in disagreement with the 
trend seen in simulations. Taking into account ellipsoidal 
collapse and the effects of large-scale environment on halo 
formation statistically (Sandvik et al. 2007) and dynami- 



cally (Desjacques 2008), the age dependent assembly bias 



can naturally arise and be more apparent for less massive 
halos. It is also noted by Desjacques (2008) that environ- 
mental density plays the determining role in the virializa- 
tion redshift of halos, and that the morphology of large- 
scale structure contributes mostly to its scatter. However, 
given a fixed halo mass, the analytical model of Desj acques] 
(2008) predicts that the median formation redshift of ha- 



los decreases with increasing environmental density, defined 
as the linear density fluctuation smoothed over 10/i~^ Mpc. 
This is inconsistent with the results from simulations. IHahnl 



et al. ( 2009 ) emphasize that for galactic-scale halos of fixed 
mass, an early formation epoch is closely related to a re- 
duced mass-growth rate at late times that is mainly due 
to the tidal effects of neighboring massive halos. As galac- 
tic halos are enhanced in filaments near massive halos, the 
dependence of their formation epoch on environmental den- 
sity is expected (see also, e.g., [Diemand et al.||2007 Wang 
et al.||2007). For cluster-scaled halos, however, Dalai et al 



(2008} show that tidal effects from large-scale structure can- 



not play an important role in the environment dependence 
of formation history. 

Gas physics allows more complicated interactions be- 
tween galaxies and their environment. Galaxies are more 
subject to ra m pressure (e.g.,|Frenk et al.| p-999) and shock 
heating (e.g., Mori fc Burkert|2000 ) in denser regions. Both 
the cooling efficiency and the gas compression rate are di- 
rectly related to environmental density. It has been ob- 
served in hydro dynamical simulations that most cold flows 
are highly anisotropic and follow dark matter filaments (e.g., 
Keres et al.||2005||Dekel fc Birnboim|2006 ). 

In this paper, we use SDSS DR7 data to examine the 
dependence of galaxy properties on their large-scale envi- 
ronment. We focus on the influence of the morphology of 
environment over and above the well-known dependence on 



environment density. This analysis may provide important 
constraints on the theories and hypotheses discussed above. 
It can also help clarify if large-scale environment should be 
included in the halo occupation distribution model (e.g., 
Blanton k Berlind |2007| [Tinker et al.|2008D . The rest of the 
paper is organized as follows. In Section 2, we describe our 
measurements of large-scale tidal field with SDSS DR7. In 
Section 3, we present our results. Comparisons with results 
from simulations with semi-analytical modeling of galaxy 
formation are shown in Section 4. Section 5 contains discus- 
sions and our conclusions. 



2 LARGE-SCALE TIDAL FIELD IN SDSS 



2.1 SDSS and NYU-VAGC 



The base data set adopted in this paper is SDSS DR7 (York 
et al.||2000| [Abazajian et al.l|2009t. S DSS used a dedicated 
wide- field 2.5 m telescope ( |Gunn et al.^^2006 ) located at 
Apache Point Observatory to obtain images in the u,g,r,i, z 
bands over an area of ~ 10,000 deg^ and spectra of ^ 1.6 
million selected objects among which galaxies amount to 
^ 1 million. The photometric data was calibrated through 
ubercalibration ( ,Padmanabhan et aH2008 ) , which used the 
overlap between adjacent imaging runs to tie together the 
photometry of all imaging observations. Fibers arranged on 
plates 1°.49 in radius were assigned with an efficient tiling 
algorithm to spectroscopic targets selected with certain pho- 
tometric criteria. Due to the fiber collision effect, only one 
spectroscopic measurement can be done for a group of tar- 
gets with separations less than 55 arsec. Overlaps of tiles can 
help to reduce the number of missed targets. The overall in- 
completeness for galaxies is ^ 6%. Our test shows that this 
incompleteness together with < 2% incompleteness caused 
by mismeasurement of redshift in spectra and bright star 
blocking is negligible for our measurements of the large-scale 
tidal field. 

Our analysis is based on the large-scale structure sam- 
ple VAGC dr72 of New York University Value Added Cata- 



logue (NYU-VAGC |Blanton et al.| ( |2005a] )). It is constructed 
from SDSS, FIRST, 2MASS, 2dFGRS and PSCz. The SDSS 
part is updated for DR7. There are ^2.5 million galaxies in 
VAGC dr72. From VAGC dr72, we build up three volume- 
limited samples in which galaxies have absolute r-band mag- 
nitude brighter than — 18 + 5 log^Lol^)' ~19 + 51og;Lo(^) ^^^ 
-20 + 51ogio(/i) respectively, {MAS, Mrl9 and Mr20 here- 
after) . We consider all spectroscopic galaxies with listed ap- 
parent magnitude brighter than r == 18. All magnitudes 
mentioned in this paper are Petrosian magnitudes, and 
are K- corrected to redshift z = 0. For Mr 18, Mr 19 and 
Mr 20, the respective number of galaxies is 26760, 53992 and 
106033. We mainly focus on the Mr 20 sample because its 
volume is large enough for accurate measurements of the 
large-scale structure geometry and it does not suffer galaxy 
bias severely since M* in the luminosity function is a little 
brighter than —20. 

The VAGC provides us information on the size and the 
stellar mass for each galaxy, so that we can directly draw 
four parameters, concentration C = R90/R50, {g — r) color, 
stellar mass M*, and the surface stellar mass density /i* = 
where R50 and R90 are defined as the radii 



M*/(27ri?io) 
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Figure 1. Left panel: The environmental density distributions for 
galaxies in eigenvalue-sign-judged halos (red solid line), filaments 
(green dashed line), sheets (blue dot-dashed line) and voids (ma- 
genta dotted line). The upper solid line shows the environmental 
density distribution for the whole sample. The smoothing scale 
for the environment of a galaxy is chosen to be the distance to 
its 3rd nearest neighbor. Right panel: Scatter plot of the environ- 
mental density 6 and ellipticity e. The solid line is the average e 
at different 6 and the dashed lines show zblcr e range around the 
average e. 




Figure 2. Contours of environmental [e, ln(l + S)] for a slice of 
SDSS DR7. The environmental density and ellipticity are encoded 
with two-dimensional color scheme shown in the insertion. The 
smoothing scale is 3h~^ Mpc. 



2.2 Tidal environment measurements 

From a theoretical point of view, the tidal environment of 
a galaxy can be characterized by the second derivatives 
of the gravitational potential, didjcj). The potential cannot, 
however, be measured directly, and needs to be calculated 
through the distribution of galaxies in redshift space. Such a 
measurement can be affected by, e.g., galaxy bias, redshift- 
space distortion and survey geometry. Here we describe our 
procedures to calculate the potential field, leaving detailed 
discussion to Section 4. 

For calculations of the potential field, we put the 
volume-limited sample into a large cubic box with side 
500/i~^ Mpc. Galaxies are assigned onto a 512^ uniform grid 
with CIC interpolation to give the number density of galax- 
ies in each cell. We remove galaxies from those cells that 
happen to cross the survey boundary, and then assign the 
average number density of galaxies to all the empty cells. 
Filling the regions outside the survey volume with the aver- 
age number density allows large voids with size comparable 
to the survey volume to be identified properly. This galaxy 
number density field is then smoothed with a Gaussian fil- 
ter with a smoothing scale Rs^ and the potential field and 
further the tidal tensor field didjCJ) at each grid point are 
calculated from the smoothed density fluctuation field. The 
tidal environment of each galaxy is obtained by linearly in- 
terpolating this discrete tidal tensor field from the mesh to 
its position. 

Without considering the specific orientation of the tidal 
tensor, its properties are fully characterized by its three 
eigenvalues Ai ^ A2 ^ A3. One way to classify the tidal en- 
vironment of a galaxy uses the signs of Ai (i = 1, 2, 3) (e.g., 
Hahn et al.|2007| [Zhang et al.|2009| ). Specifically, one can de- 
fine a point as (i) "cluster" , if all three eigenvalues are posi- 
tive; (ii) "filament" , if one eigenvalue is negative; (ii) "sheet" , 
if two eigenvalues are negative; and (iv)"void", if all the 
three eigenvalues are negative. These definitions can show 
clearly the filament-node nature of large-scale structure at 
least in a two-dimensional cut. However, this classification is 
highly degenerate with a classification by density. As shown 
in the left panel of Figure 1, the overdensity distributions 
of points in these different categories are very different, and 
thus the dependence of galaxy properties on such "tidal en- 
vironments" may mainly reflect a density dependence. To 
address the dependence of galaxy properties on tidal prop- 
erties other than overdensity, we here describe the tidal en- 
vironment of galaxies with parameters describing ellipticity 
e and prolateness p. Specifically, we define e and p as follows 



Ai — A3 

3 + ^ ' 



Ai + A3 — 2A2 



containing 50% and 90% of the Petrosian flux, respectively. 
In VAGC, calculations for the stellar mass M* are based on 



five-band photometric data and redshift. Li & White (2009) 



compare this stellar mass measurement with spectroscop- 
ically based stellar mass measurements by [Kauffmann et| 



al. (2003), and show that the two results agree very well. 



We also add the amplitude of the 4000 A break Dn(4000) 
into our analysis as a clear indicator of stellar population 
age, from the MPA/JHU SDSS data {littp://www.mpa- 
garching.mpg.de/SDSS/). There are 105034 galaxies in 
Mr20 that have counterparts in MPA/JHU data set. 



where ^ = Ai + A2 + A3. We find that these definitions, with 
the number 3 added to the denominators, can minimize the 
correlations between e (p) and 6, as seen in the right panel 
of Figure 1. Figure 2 shows a 2-dimensional color rendering 
of [e, ln(l + 6)] in a 2-dimension cut through the Mr20 sam- 
ple. One can see that the cluster-filament-cluster structure 
is very well represented. High density and high ellipticity 
bridges link high density and low ellipticity clusters. In the 
following analysis, we mainly focus on the e-dependence of 
galaxy properties, since the ^-dependence has been well ex- 
plored in previous work. 
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Figure 3. Galaxy color contours in the ellipticity-density plane for the SDSS DR7. The fixed smoothing scales are \h~^ Mpc (left 
panel), 2h~^ Mpc (middle panel) and S/i"-*^ Mpc (right panel), respectively. The interval of the color contours is taken to be O.OScr. 
The black rings are the uncertainty contours of galaxy color with values O.OScr, O.lcr, O.lScr and 0.2cr from inside out, where a = 0.135 
is the standard deviation in galaxy color for the whole Mr 20 sample. The horizontal black lines are plotted so that the behavior of 
correlation/anticorrelation with ellipticity can be seen. 
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Figure 4. Galaxy color contours similar to Figure 3. The smoothing scales are lh~^ Mpc (left panel), the 'adaptive smoothing scale' 
(middle panel) and ?th~^ Mpc (right panel), respectively. Here the 'adaptive smoothing scale' scale applied to smooth the environment 
of a galaxy is defined as the distance to its 3rd nearest neighbor. The black contours shown in the left and right panels are the contours 
of ln(l + (5o) with 5o the environmental density measured by this same optimal adaptive smoothing. 



3 THE TIDAL DEPENDENCE OF GALAXY 
PROPERTIES 

3.1 Tidal dependence of galaxy color 

As can be seen from Figure 2, giant filaments link clusters 
together. These filaments have similar densities to the out- 
skirts of clusters. By introducing the additional parameter 
e, however, these two regions can be separated clearly. 

The main purpose of our study is to investigate if galaxy 
properties depend on the morphology of their large-scale 
environment in addition to the known density dependence. 
To do that, we consider galaxy properties in the [e, ln(l + ^)] 
plane, and analyze the correlations of these properties with 
e and ln(l + ^). Specifically, we divide the [e,ln(l + ^)] plane 
into 128^ cells, where e G [0, 1.2] and the typical range of 
ln(l + (5) is (—1, 4) with its precise range varying somewhat 
from one smoothing scale to another. We then extend the 
plane to a periodic 256^ grid by reflecting the original cells 



with respect to e = and a line of the lowest value of ln(l + 
5). At each cell, we calculate the total number of galaxies 
A/'[e, ln(l+^)] and the sum of the considered galaxy property 
C[e, ln(l + 5)]. Smoothing is then applied to both C and N 
with a Gaussian smoothing function Cb, where subscript 5 
denotes that the smoothing scale is 5 cells, and the average 
{C)s is obtained at each grid point by {C)s = Cs/Ns, where 
the subscript s represents the smoothed quantities. 

To quantify the statistical significance of the e- 
dependence of galaxy properties, we need to estimate the 
statistical errors for the smoothed quantities properly. Our 
sample contains ^0.1 million galaxies. When they are dis- 
tributed onto 128^ cells in the [e,ln(l + S)] plane, we ex- 
pect that Poisson fluctuations due to the limited number 
of galaxies associated with each cell will dominate the sta- 
tistical errors. We thus estimate the error in each cell as 
crV^A^/A/"!, where a is the standard deviation of the consid- 
ered galaxy property over the whole sample, and A^i and N2 
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Table 1. Average value and standard deviation of different galaxy 
properties for Mr 20 



Property 


(g-r) 


-Mr 


Dn4000 


/i* 


R90/R50 


M* 


mean 
a 


0.663 
0.135 


20.63 
0.481 


1.676 
0.259 


9.076 
0.376 


2.706 
0.433 


10.30 
0.256 



are the smoothed cell count with the smoothing functions 
G5 and Gi, respectively. In Table 1, we list the values of a 
for different galaxy properties for Mr 20. 

Figure 3 shows the results for the color (g — r), where 
a — 0.135 as derived from the sample Mr 20 and shown in 
Table 1. The three different panels correspond to three dif- 
ferent smoothing scales to represent the large-scale environ- 
ment, Rs = lh~^ Mpc, 2h~^ Mpc and 3h~^ Mpc from left 
to right. The (g — r) values are shown in color in each panel 
with the scale shown to the right of Figure 3. For the super- 
posed white contour lines, the interval between the adjacent 
lines is 0.05cr. The inside-out black lines show the error con- 
tours of 0.05cr, O.lcr, 0.2(7 and 0.3cr, respectively. 

The dependence of (g—r) on the density S is clearly seen 
in Figure 3. The slope of the (g — r) contours reflects whether 
the ellipticity e also affects the galaxy color in addition to S. 
We can see a weak correlation between (g — r) and e, and this 
correlation is scale-dependent. For Rs — lh~^ Mpc, an anti- 
correlation is seen. This reverses with increasing Rs^ and is 
positive for Rg — 3h~^ Mpc and larger. For Rs = 2h~^ Mpc, 
the correlation is nearly null. 



3.2 Scale-dependence of environmental effects 

From Figure 3, we see that galaxy colors depend strongly 
on environmental density 6 and weakly on environmental 
ellipticity e. This behavior are scale-dependent. The ques- 
tion arises then what smoothing scale one should adopt in 
order to reveal the most fundamental dependence of galaxy 
properties on their environment. 

Previous studies have shown that galaxy properties are 
inffuenced by their environmental densities mainly locally 



jKauffmann et al.|[2004l [Blanton k Berlind ||2007| |Park"et 
|al.||2007l). On the other hand, if we simply smooth the en- 
vironments on a scale larger than expected for any physical 
impact, we may still see an apparent environment- density 
dependence which arises merely due to the correlation be- 
tween the environmental density on the considered smooth- 
ing scale and that on the scale where environmental effects 
are the strongest. This can also apply to the apparent e- 
dependence of galaxy properties shown in Figure 3. 

The scale-dependent switch from anti-correlation to 
positive-correlation with e for galaxy colors seen in Fig- 
ure 3 indicates to us that the best smoothing scale should 
be around the transition scale where no e-correlations are 
observed. Furthermore, this transition scale should corre- 
spond to the scale where the density dependence of galaxy 
properties is the strongest. In Appendix A, we discuss fur- 
ther this scale-dependence of the environmental-density in- 
ffuence on galaxy properties. We find that the scale where 
the environmental-density dependence of galaxy proper- 
ties is the strongest is somewhat different for high and 
low density regions with the values of ^ 1.5/i~^ Mpc and 



2.5h~^ Mpc, respectively. This difference suggests that an 
adaptive smoothing scale may be more appropriate. A nat- 
ural choice is the characteristic distance between the con- 
sidered galaxy and its neighbors. Our analysis show that 
for Mr 20 sample, the distance to a galaxy's 3r(i- nearest- 
neighbor is an optimal adaptive smoothing scale, which has 
a distribution peaked at 2h~^ Mpc and is close to the peak 
scale where galaxy properties and environment density have 
the strongest correlation in both high and low density re- 
gions. For Mr 19 and Mr 18 samples, this optimal scale should 
correspond to the distance to a galaxy's 6t/i-nearest- neigh- 
bor and 8t/i-nearest-neighbor respectively, due to the fact 
that the galaxy density is higher. 

In Figure 4, we show {g — r) contours as in Figure 3. 
The right and left panels are for fixed smoothing scales of 
lh~^ Mpc and 3h~^ Mpc, respectively. The middle panel 
shows results using the 3rd-nearest-neighbor distance as an 
adaptive smoothing scale. In comparison with the middle 
panel of Figure 3, the null-dependence on e is more cleanly 
seen for the adaptive smoothing. To investigate if the de- 
pendence on e seen in the left and right panels of Figure 
4 (the same as those of Figure 3) is physical or due to the 
correlations between (e, 6) at the considered scale and 6 at 
the 'optimal adaptive smoothing scale', we calculate the lat- 
ter correlations. Specifically, we define (e, S) at the consid- 
ered smoothing scale and at the optimal adaptive smooth- 
ing scale as (ca^Sa) and {eo,So), respectively. We compute 
the ^o-contours in the (ca^Sa) plane, which are shown as 
black lines in the left and right panels of Figure 4. It is 
seen that the ^o-contours align with the (g — r)-contours 
(white lines) very well. This demonstrates that the observed 
e-dependence is mainly attributable to the correlation be- 
tween (ca^Sa) and So- Some deviation between the (g — r)- 
contours and the ^o-contours are seen in high density re- 
gions. These are, however, insignificant compared to the sta- 
tistical errors. Together with the null-correlation with e for 
the optimal adaptive smoothing, we conclude that no signif- 
icant physical dependence of galaxy colors on e is detected 
by our analysis. 

3.3 Other galaxy properties 

We have studied the dependence of galaxy color on (e, ^) 
carefully. Here we examine the tidal dependence of other 
galaxy properties. 

Galaxy properties can be divided into three classes: (1) 
star-formation-related properties such as galaxy color and 
Dn(4000); (2) morphology-related properties such as concen- 
tration; and (3) extensive quantities such as the total stellar 
mass M* and galaxy size. Previous studies have shown that 
star-formation-related properties depend directly on the en- 
vironmental density (Kauffmann et al^|2004 Blanton et al 



2005b Christlein &: Zabludoff 2005), but morphology- 



related ones are not independently correlated with environ- 
mental den sity ( [Park et al.||2007| [Ball et al.||2008| [Bamford 
eraL][2 008l. 

Here we analyze the (e, ^)-dependence of different 
galaxy properties. Figure 5 presents the environmental de- 
pendence of (M*,Dn4000). The left panel shows a scatter 
plot of (M*,Dn4000). The red lines are contours of con- 
strained mean value of ln(l + 5o)- From the scatter plot, 
we see that Dn4000, which reflects stellar population ages 
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Figure 5. Left panel: The scatter plot of (M*, Dn(4000)). The 
ln(l -\- 5o) contours are shown in red. Right panel: The contours 
of Dn(4000)-{Dn(4000)>, where {Dn(4000)> is the average value 
of Dn(4000) given stellar mass. 
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Figure 6. Left panel: The scatter plot of (Dn(4000), C). The 
ln(l + 5o) contours are shown in red. Right panel: Contours of 
C — (C), where {C) is the average value of C given Dn(4000). 



and therefore star formation histories, correlates with total 
stellar mass M* . The directions of the red contours indicate 
that both depend on environmental density as measured by 
our optimal adaptive smoothing. In the right panel of Fig- 
ure 5, we show the (e, 5) dependence for Dn4000 given the 
stellar mass M* . Specifically, the contours of the deviation 
Dn4000 - (Dn4000) are shown, where (Dn4000) is the aver- 
age value of Dn4000 in stellar mass bins. The environmental- 
density dependence of the contours are clearly seen, which, 
in agreement with the trend seen in the left panel, demon- 
strates the extra density dependence of 0^4000 in addition 
its dependence of M*. On the other hand, no dependence on 
the shape parameter e is detected. 

In Figure 6, we show a similar analysis for the con- 
centration parameter C = Rqq/Rbq. The scatter plot of 
(Dn4000, C) is shown in the left panel with the ln(l + 6o) 
contours superposed on it. Two clumps are seen in the scat- 
ter plot with one at low value of Dn4000 and small value of 
C, and the other at high 0^4000 and large C. The lower- 
left clump corresponds to galaxies with recent star forma- 
tion, and the upper-right clump is associated with galaxies 
with quenched star formation. The ln(l + ^o) contour lines 
are nearly vertical, showing no detect abe independent cor- 
relation between C and the environmental density. In the 
right panel of Figure 6, we show the C — (C) contours in 
the [eo,ln(l + So)] plane, where (C) is the average value of 
C given Dn4000. The chaotic pattern seen here reflects the 
lack of any environmental dependence of the concentration 
C other than its So dependence through the correlations with 
Dn4000. 

We also analyze the environmental dependence of 
galaxy size defined as i?5o*d in unit of h~^kpc, where d 
denotes the distance to a galaxy. The results are shown in 
Figure 7, again using our optimal adaptive smoothing. It 
is seen that in high density regions, galaxy size has no de- 
tectable correlation with either environmental density or el- 
lipticity. In low density regions, on the other hand, there is 
a notable correlation between galaxy size and Co in addition 
to its correlation with So. We will show in section 4.3 that 
in low density regions, the survey boundary effects can bias 
the So and Co measurements to lower values. This can in- 
duce an artificial correlation with Co, which is considerable 
for small- volume samples such as Mr 18. For Mr 20, how- 
ever, the survey volume is large and the boundary effects 
are insignificant. Thus the correlations seen in the lower left 




3.5 



3.15 



Figure 7. Contours of galaxy size defined as R^o * Rz in unit of 
for SDSS. The boundary is the 0.2cr uncertainty contour, where 
a is the standard deviation of galaxy size for the whole sample. 
The optimal adaptive smoothing is applied. 



region of Figure 7 may indicate that both the density and 
the ellipticity of the environment, can affect the galaxy size. 
Galaxies are systematically smaller in lower (So, Co) regions. 
Due to the still limited statistics of SDSS, the environmental 
Co dependence of galaxy size needs to be investigated further 
with future large surveys that can provide much improved 
statistics. 

Our analysis so far has measured ellipticity and density 
on the same scale. The optimal adaptive smoothing is chosen 
to be the distance from a galaxy to its 3rd-nearest neighbor. 
This scale has a distribution peaked between 1.5h~^Mpc and 
2.5h~^Mpc, depending on environmental density. In order to 
show the influence of larger scales, especially the influence of 
the large filaments linking cluster as seen in Figure 2, we also 
analyze the environmental dependence of galaxies properties 
on ellipticity and density measured on two different scales. 
We fixs the smoothing for density to the optimal adaptive 
smoothing and we change the smoothing scale for ellipticity. 
For no smoothing scale larger than the 2h~^Mpc do we see a 
correlation of galaxy properties with environment ellipticity. 
We thus conclude that large scale structure beyond galaxy 
groups has no additional influence on individual galaxies. 
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We have so far focused our analysis on the eUipticity 
e of large-scale environment. However, we have also looked 
for the dependence of galaxy properties on the environmen- 
tal prolateness p. We find that smoothing according to the 
distance to the 3rc/-nearest-neighbor of a galaxy, all of the 
considered galaxy properties are independent of p. 



4 COMPARISON WITH A SEMI-ANALYTIC 
MODEL 

In semi-analytic models of galaxy formation (SAM), our best 
understanding on dark matter halo formation is combined 
with the complex baryonic physics associated with galaxy 
formation as represented by physically/observationally mo- 
tivated prescriptions involving a set of parameters. While 
earlier SAMs were based on using the extended Press- 
Schechter theory to construct dark matter halo merger trees, 
modern SAMs take full advantage of high resolution N-body 
simulations to follow the formation and evolution of dark 
matter halos. Galaxies are then formed within these halos 
by modeling physical processes such as cooling, condensa- 
tion, star formation, supernova and AGN feedback. Here we 

Blaizot] 



analyze the SAM galaxy catalog from |De Lucia 



I (2007), for which the dark matter distribution was taken 
from the Millennium Simulation (Springel et al. 2005). We 
compare the results with those from the SDSS. To test obser- 
vational effects, we also construct mock SDSS DR7 catalogs 
from the SAM data ( Li et al.|2007 ) . Specifically, we address 
the following issues 

(i) a comparison of the environmental-density depen- 
dence of galaxy properties in the SAM and the SDSS; 

(ii) a comparison of the behavior of the tidal dependence 
of galaxy properties; 

(iii) the influence of various observational effects, such as 
redshift-space distortion and survey geometry on the results 
of our tidal analysis. 



4.1 Scale dependence of the color-density relation 
in the SAM 

To analyze the scale dependence of the color-density rela- 
tion, we calculate the correlation coefficient between galaxy 
color (g—r) and environmental density S at different smooth- 
ing scales. The correlation coefficient is defined as 

(A(3 - r)ASo} 



where A{g — r) = {g — r) — {{g — r)) and A6o = So — (So). 
The results are shown in Figure 8. We divide galaxies into 
two groups of high (solid lines) and low densities (dashed 
lines) as described in Appendix A (the green lines in Figure 
Al). The red, blue and green lines are, respectively, for the 
SDSS, the full volume of the SAM, and for the mock SDSS. 
First, we see similar trends for the scale dependence 
of the color-density correlation in the SAM and the SDSS 
data. The SAM 'optimal smoothing scale' for the strongest 
correlation is also consistent with that of SDSS, which is 
~ lh~^ Mpc for high density regions and ^ 2h~^ Mpc for 
low density regions. On the other hand, the amplitude of 
the correlation coefficient is significantly larger for the SAM 



SDSS high density 
SDSS low density 
mock low density 
mock high density 
SAM low density 
SAM high density 




2 3 4 5 

Smoothing Scale (Mpc h"^) 



Figure 8. The color-density correlation coefficient as a function 
of smoothing scale. The red, blue and green lines give results 
for SDSS, SAM, and mock catalog, respectively. The solid and 
dashed lines are the results for galaxies in high and low density 
regions, respectively (see Appendix A). The filled circles of dif- 
ferent color represent the results for optimal adaptive smoothing 
of the corresponding samples. 



than for the SDSS. The SAM predicts an environment- 
density dependence of galaxy colors stronger than is ob- 



served, in accordance with other studies (e.g., Springel et 



ar||2005| ICbil et al.||2008| [Font et al.||2008| |Cowanfc Ivezic 
2008||Guo et aL|2011t . We note that while the tidal interac- 
tions of dark matter halos are naturally accounted for in the 
SAM, tidal disruption of galaxies is not. For massive clus- 
ters, however, it is expected that their strong tidal forces 
may destroy galaxies (e.g., Henriques et al. 2008). As those 
satellite galaxies are mostly red, such disruptions can de- 
crease the level of clustering of red galaxies, and therefore 
reduce the discrepancies in color-density relation between 
the SAM and the observations in high density regions. Fur- 
ther, this can also alleviate the overly- fast-growth problem 
for central galaxies in the SAM ( Wang fc Jing|2010 ). It has 
also been argued that hot gas in satellite galaxies may not 
be immediately stripped away by ram pressure when they 
merge into a larger system as often assumed in SAMs. There- 
fore a reservoir of gas can be replenished to maintain the star 
formation for a relatively long period of time (e.g., [Fbnt et 



al][2008l [Coil et al.|[2008l [Weinmann eraL][2006l |Guo et al. 



2011 ). In low density regions, the SAM also predicts stronger 



color-density correlations than seen in the SDSS, which re- 
flects the smaller color spread in the SAM compared to the 
SDSS. 

The green lines in Figure 8 are the results from the 
SDSS DR7 mock catalogue constructed from SAM data. The 
differences between the green and blue lines reflect observa- 
tioal effects which will be discussed in ^4.3. 
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Figure 9. Galaxy color contours similar to Figure 4, but for the SAM. Here the black contours are the color contours for mock data. 
The smoothing scales are lh~^ Mpc (left panel), distance to the Srd-nearest-neighbor (middle panel) and ^h~^ Mpc (right panel), 
respectively. The outer boundary is the uncertainty contour of galaxy color with the value of O.lcr, where cr = 0.157 is the standard 
deviation of the galaxy color distribution for the whole SAM sample. 



4.2 Tidal dependence of galaxy properties in 

SAM 

To study the tidal dependence of galaxy properties, we per- 
form the same analysis for the SAM data as shown in §3.2 
for the SDSS. The results are shown in Figure 9. The black 
lines here represent results from the mock catalogue, and 
will be discussed in §4.3. The three panels correspond to 
the smoothing with Rs — lh~^ Mpc, with the 3r(i- nearest- 
neighbor distance, and with the Rs — 3h~^ Mpc. It is 
seen that the contour behaviors are qualitatively similar to 
those of Figure 4 for the SDSS. In the SAM, the distance 
to the 3r(i-nearest neighbor also corresponds to the 'peak 
scale' of the color-density dependence, at which, the tidal 
e-dependence of galaxy properties is nearly null. The den- 
sity So at the optimal adaptive smoothing scale plays the 
dominant role for the environmental effects. The weak nega- 
tive/positive {g — r)-e correlations at smaller/larger smooth- 
ing scales again reflect correlations between (e, 6) at the 
smoothing scale and So. 

A notable difference between Figure 9 and Figure 4 
is the density dependence of color, which is significantly 
stronger in Figure 9 than that shown in Figure 4. This is 
consistent with the result shown in Figure 8, which also show 
that the SAM predicts a considerable stronger color-density 
correlation than seen in observations. 



4.3 Observational effects 

To make a realistic comparison between the results for the 
SAM and the SDSS, different observational effects must be 
taken into account. We thus build mock catalogs from the 



SAM data following Li et al. (2007). We first test the influ- 



ence of fiber collisions which reduce the number of redshift 
measurements for galaxies in close pairs by about 6%. Be- 
cause our smoothing scale is nearly always larger than 55^^ 
fiber collisions do not affect our analysis significantly. 

The most notable influences from obervations are the 
redshift distortion effect and the effect from the survey ge- 
ometry. The redshift-space distortion affects mainly the high 



density regions and plays a role of "mixing". In redshift 
space, galaxies in massive clusters present a Finger-of-God 
configuration along the line of sight. On the other hand, 
galaxies in filaments surrounding a massive cluster have a 
tendency to be flowing into the cluster, which squashes the 
filaments along the line of sight. In other words, galaxies 
in filaments and galaxies in massive clusters are somewhat 
mixed in redshift space. Thus the color scatter in a given 
density bin becomes larger in redshift space, resulting in 
a lower color-density correlation coefficient in high density 
regions, as shown by the green solid line in Figure 8. Fur- 
thermore, the redshift-space distortion makes the measured 
environmental density systematically lower than that mea- 
sured in real space. The ellipticity calculated in the red- 
shift space is systematically lower/higher for regions with 
high/low e in real space. That is, filaments become rounder 
and clusters tend to be more anisotropic in redshift space. 
Such effects cause color mixing in the [e,ln(l + S)] plane, 
and thus lead to flatter color contours with respect to the 
ellipticity e and weaker correlations with 6 (see the black 
contour lines in Figure 9). 

In low density regions, measurements of environmental 
density and ellipticity are mainly affected by survey geome- 
try. Both are biased to lower values. In the left panel of Fig- 
ure 10, we show the value of e with respect to the distance to 
the survey boundary for SDSS sample Mr 20. Different lines 
represent different environmental density ranges with longer 
dashes corresponding to higher densities. Physically, we do 
not expect any correlations between e value and the dis- 
tance to the survey boundary. For high-density lines, they 
are indeed nearly flat. For low- density lines, however, the 
measured e is systematically lower for galaxies closer to the 
survey boundary. The effect is significant for galaxies that 
are within 20h~^ Mpc of the boundary. This is clearly an ob- 
servational effect due to padding the regions outside the sur- 
vey boundary with the average number density. Such a mea- 
surement bias on e can lead to artificial correlations between 
galaxy properties and environmental ellipticity. The signif- 
icance of the effect depends on the the fraction of galaxies 
that are within 20h~^ Mpc of the boundary. Volume- limited 
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samples with smaller volume are affected more. In the right 
panel of Figure 10, we show the ratio of color in high e re- 
gions to the color in low e regions vs. minus the absolute 
magnitude of galaxies -Mr. Here high and low e are defined 
as above the la level and below the —la level with respect to 
the average e (see the right panel of Figure 1). We present 
results for different SDSS volume-limited samples, Mr 18, 
Mr 19 and Mr 20 from smaller to larger volume. For Mr 18, 
we see the decrease of the color ratio with respect to -Mr 
of galaxies. That is, for relatively faint galaxies, the ones in 
low e regions tend to be bluer than those in high e regions. 
For bright galaxies, the ratio is close to 1. For Mr 19 and 
Mr 20, there is no detectable dependence of the color ratio 
on -Mr, and the ratio is always close to 1. To test if the Mr- 
dependence of the color ratio seen in Mr 18 is physical or due 
to the boundary effect, we construct subsamples of galaxies 
with distances less than 110h~^ Mpc from Mrl9 and Mr20, 
respectively. This cut corresponds roughly to the distance 
limit of Mrl8. The results for the Mrl9 and Mr20 subsam- 
ples are shown by unfilled green triangles and unfilled red 
squares, respectively. The apparent correlation with Mr of 
galaxies is then seen. This demonstrates that the correlation 
seen for Mr 18 sample is largely due to the boundary effect. 
As shown in the left panel of Figure 10, in low density re- 
gions, the e value is biased to a lower value for galaxies closer 
to the survey boundary. Because galaxies in low density re- 
gions tend to be bluer than those in high density regions, 
this bias leads to an increase in the number of relatively blue 
galaxies in low e regions, and thus increases the color ratio 
between high e and low e regions. The increase of the color 
ratio is more significant for fainter galaxies because the ob- 
served ones reside in a smaller volume around the observer 
and thus suffer more from the boundary effect. 

In Figure 11, we show color contours in the [eo,ln(l + 
6o)] plane for Mr 18. For Mr 18, the optimal adaptive smooth- 
ing scale for a galaxy corresponds to the distance to its 8th 
nearest neighbor. Comparing to the middle panel of Figure 
4, we see a notable correlation between color and ellipticity 
in low-density regions. This however, is largely due to the 
survey boundary effect that biases both the environmental 
density and ellipticity toward lower values in low density 
regions. 



5 CONCLUSION AND DISCUSSION 

In this paper we examine the tidal dependence of galaxy 
properties in the NYU-VAGC sample and compare it with 
the predictions of the SAM. To separate the environment- 
morphology dependence from the environmental-density de- 
pendence, we construct the density fluctuation 6, the elliptic- 
ity e and the prolateness p from the three eigenvalues of the 
tidal tensor of the potential field calculated from the spatial 
distribution of galaxies, where e and p are nearly indepen- 
dent of ^. It should be noted that this potential field can be 
biased with respect to the potential field from the underlying 
dark matter distribution. The well-known galaxy bias refers 
to the ratio of the density fluctuation amplitude of galaxy 
distribution and that of the dark matter distribution. 

For Mr 20, the volume-limited sample with absolute 
magnitude brighter than —20, our results show that while 
the environmental-density dependence of galaxy properties 
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Figure 10. Left panel: The average environment ellipticity for 
galaxies as a function of their distance to the survey boundary. 
The optimal adaptive smoothing has been applied. Different lines 
represent the results for different bins of log{l + (5o) ranging from 
log(l + (5o) = —2.0 to ln(l + 5o) = 4.0. Longer-dashed lines repre- 
sent the results for higher density bins. Right panel: The ratio of 
the mean (g — r) in high ellipticity regions to that in low ellipticity 
regions for Mr 20 (red filled circle), Mr 19 (green filled circle) and 
Mr 18 samples (blue filled circle), respectively. The high and low 
ellipticity regions are defined as the regions above the +lcr dashed 
line and below the —la dashed line shown in the right pane of 
Figure 1. The results for subsamples of Mr 20 and of Mr 19 cut 
at 110h~^ Mpc are shown by red open squares and green open 
triangles, respectively. 
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Figure 11. Galaxy color contours for SDSS Mr 18. The environ- 
mental smoothing scale for a galaxy is taken to be the distance 
to its 8t/i-nearest neighbor. The boundary is the 0.3cr uncertainty 
contour, where a = 0.168 is the standard deviation of galaxy color 
of the whole Mr 18 sample. 



is indisputable, the dependence on e and p is rather weak, 
except perhaps for galaxy size. It is further shown that such 
weak correlations depend on the smoothing scale consid- 
ered for the environment. There exists a particular scale, 
^ 2h~^ Mpc, at which, the environmental- e (p) dependence 
nearly vanishes. This corresponds to the optimal adaptive 
smoothing scale where the environmental-density depen- 
dence of galaxy properties is the strongest. It varies from 
high density regions to low density regions, and corresponds 
well to the distance to each galaxy's Srd-nearest-neighbor 
for Mr 20 sample. At smoothing scales larger /smaller than 
this 'optimal smoothing scale', the weak correlations be- 
tween galaxy properties and e are negative/positive. We 
demonstrate, however, that these correlations result mainly 
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from correlations between (e, ^) at the considered smoothing 
scale and the environmental density 5o at the optimal adap- 
tive smoothing scale. In other words, we find no physical 
influence of environment morphology on galaxy properties. 
Furthermore, it is 5o that plays the dominant role for the 
environmental effects, and the apparent density dependence 
on other smoothing scales is largely due to the correlations 
between 5 and 5o. This indicates that, galaxy properties in 
Mr 20 are affected mainly by their nearby environments. Our 
analysis shows that galaxy size is independent of environ- 
ment in high-density regions. In low-density regions, some 
correlation with Co in addition to the dependence on 5o is de- 
tected. Galaxies in lower (^o, Go) regions tend to be smaller. 
Such correlations need to be further explored with future 
observations. 

With NYU-VAGC M^IQ, [Blanton fc Berlind | ( [2007 ) 
analyze galaxies in groups, and find that galaxy proper- 
ties depend mainly on the properties of their host groups, 
and are not affected independently by the environment on 
larger scales(> lh~^ Mpc). As they only consider galaxies in 
groups, it is appropriate to compare their results with ours 
in high density regions. From Figure 8, we can see that the 
'peak scale' in high density regions is ^ lh~^ Mpc, consis- 
tent with the typical group scale in their studies. Therefore 
our results are in agreement with theirs concerning the dom- 
inant importance of the nearby environment. This in turn 
validates the current halo occupation distribution model, 
which assumes that galaxy properties depend only on those 
of their parent halos with no additional influence from larger 
scale. Lee k, Lee] ( ,2008J also examines the tidal dependence 
of galaxy properties. They found that there is a tendency 
for elliptical galaxies to be preferentially situated in low 
e regions with density in the range 0.5 ^ 5 ^ 1.06. For 
—0.3 ^ 5 ^ 0.1, spiral galaxies are likely in regions with 
large e. We note that their considered scale for environmen- 
tal effects is ^ 400/64 r^ 6h~^ Mpc. Thus their results can 
loosely be compared with our analysis for Rs = 3h~^ Mpc 
(but note the different 6 range). With the correspondence 
between color and morphology of galaxies, the positive cor- 
relation of galaxy color with e seen in right panel of Figure 
4 is qualitatively consistent with the trend found byl Lee fc] 



Lee 



(2008) although their definition of e is different from 



ours. On the other hand, our analysis indicates that this e- 
dependence is largely due to the correlations between e and 
6o and therefore has little physical significance. 

Comparing the SDSS results with those from our SAM, 
we see qualitative agreement between the two, although 
the SAM predicts a stronger dependence of galaxy prop- 
erties on environmental density. Our results are also in ac- 



cordance with the theoretical study of Desjacques (2008), 



which showed that environmental effects on dark matter 
halo formation mainly reflect environmental density, and are 



influenced little by e and p. Hahn et al. (2009) found that 



the anti-correlation between the formation epoch of galactic 
halos and their environment density is, in fact, mainly at- 
tributable to tidal suppression by neighboring massive halos. 
Because of the enhanced environment density near massive 
halos, this suppression shows up as an assembly bias, and 
provides a possible explanation for galaxy variation with en- 
vironment density. On the other hand, jHahn et al.| ( |20d9| 
chose the smallest eigenvalue A3 of the tidal tensor as an 
indicator of the strength of tidal field, and they found a 



stronger dependence of halo formation epoch on A3 than 
on environment density. Our results manifest a strong de- 
generacy between A3 and environment density. Just as for 
ellipticity, A3 can be geometrically more strongly correlated 
with environment density on the optimal smoothing scale 
than on a larger scale. Further work is needed to clarify 
whether the dependence of halo formation history on en- 
vironment morphology can be attributed to this geometry 
effect, or whether baryon physics actually make difference 
between halos and galaxies. 

By constructing mock catalogs from SAM data, we have 
examined the infiuence of observational effects on our anal- 
ysis. Redshift-space distortion is significant in high density 
regions. It mixes galaxies in filaments and galaxies in clus- 
ters causing an increase of color scatter at given density, 
and therefore reduces the color-density correlation coeffi- 
cient considerably. The survey boundary affects environment 
measurements mainly in low density regions, biasing both 
the environmental density and the ellipticity to lower val- 
ues. The effect is notable for galaxies within 20h~^ Mpc of 
the survey boundary. While it is negligible for Mr 20, the 
effect is considerable for Mr 18 which occupies a relatively 
small volume and therefore has a large fraction of its galax- 
ies close to the boundary. It induces an artificial correlation 
between galaxy properties and the environmental ellipticity 
for Mr 18. 

In summary, our analysis of SDSS data shows that 
in addition to environment density, there is no significant 
further dependence of galaxy properties on the morphol- 
ogy of large scale structure. Geometrically, both ellipticity 
and density on one smoothing scale correlate strongly with 
each other, and with ellipticity/density on other smooth- 
ing scales. If the smoothing scale is not chosen properly, 
an appeared dependence of galaxy properties on both envi- 
ronment ellipticity and environment density arises which is 
merely due to geometry. We find that for the optimal adap- 
tive smoothing scale, the dependence on density is maxi- 
mized and the dependence on ellipticity and prolateness is 
null. 



ACKNOWLEDGMENTS 

We thank Cristiano Porciani and Guinevere Kauffmann for 
detailed discussions in the early stages of this project, and 
Cheng Li for his help on dealing with the SDSS data. This 
research is supported in part by the NSFC of China un- 
der grants 10773001 and 11033005, 11173001 and the 973 
program No. 2007CB815401, as wefl as by the ERG under 
Advanced Grant 246797 "Galformod" . 



REFERENCES 

Abazajian K., Adelman-Mccarthy J.K., Agueros, M., 2009, 

ApJS, 182, 543 
Bah N.M., Loveday J., Brunner R.J., 2008, MNRAS, 383, 

907 
Bond J.R., Myers S.T., 1996, ApJS, 103, 1 
Bamford S.P., Rojas, A.L., Nichol R.C., 2008, MNRAS, 

391, 607 



© 2010 RAS, MNRAS 000,[T}{TT] 



On the tidal dependence of galaxy properties 11 



Blanton M.R., Schlegel D.J., Strauss M.A., 2005, ApJ, 128, 

2562 
Blanton M.R., Eisenstein D., Hogg D.W., 2005, ApJ, 629, 

143 
Blanton M.R., Berlind A. A., 2007, ApJ, 664, 791 
Bond J.R., Efstathiou G., Kaiser N., 1991, ApJ, 379, 440 
Christlein D., Zabludoff A., 2005, ApJ, 621, 201 
Coil A., Newman J.A. & Croton D., 2008, ApJ, 672, 153 
Cowan N.B., IveziC Z., 2008, ApJ, 674, 13 
Dalai N., White M., Bond J. R., Shirokov A., 2008, ApJ, 

687,12 
Dekel A. Birnboim Y., 2006, MNRAS, 368,2 
De Lucia C, Blaizot J., 2007, MNRAS, 375, 2 
Desjacques V., 2008, MNRAS, 388, 638 
Diemand J., 2007, ApJ, 667, 859 
Dressier, A., 1978, ApJ, 223,765 
Croton D.J., Gao L., White S.D.M., 2007, MNRAS, 374, 

1303 
Gao L., Springel V., White S.D.M., MNRAS, 363, 66 
Gao L., White S.D.M., 2007, MNRAS, 377, 5 
Gunn J.E., Siegmund W.A., Mannery E.J., 2006, AJ, 131, 

2332 
Guo, Q., White, S.D.M., Boylan-Kolchin, M., & De Luica, 

G., 2011, MNRAS, 413, 101 
Frenk C.S., White S.D.M., Bode P., ApJ, 525, 554 
Font A.S., Bower R.G., McCarthy I.G., MNRAS, 389, 1619 
Hahn O., Corollo M., Porciani C, 2007, MNRAS, 381,41 
Hahn O., Porciani C, Dekel A., 2009, MNRAS, 398, 1742 
Henriques B. M., Bertone, S., Thomas P.A., 2008, MNRAS, 

383, 1649 
Hubble E.P., 1936, The Realm of the Nebulae. New Haven, 

CT: Yale Univ. Press 
Harker G., Cole S., Helly J., 2006, MNRAS, 367, 1039 
Jing Y.P., Suto Y., Mo H.J., 2007, ApJ, 657, 665 
Kauffmann G., Heckman T.M., White S.D.M. 
Kauffmann G., White S.D.M., Heckman T. M., 2004, MN- 
RAS, 353, 713 
Keres D., Katz N., Weinberg D.H., 2005, MNRAS, 363, 2 
Keselman J.A., Nusser A., 2007, MNRAS, 382, 1853 
Kim H.S., Baugh CM., Cole S., 2009, MNRAS, 400, 1527 
Lacey C, Cole S., 1993, MNRAS, 262, 627 
Lee J., Lee B., 2008, ApJ, 688, 78 

Li C, Kauffmann G., Jing Y.P, 2006, MNRAS, 368, 21 
Li C, Jing Y.P, Kauffmann G., 2007, MNRAS, 376, 984 
Li C, White S.D.M., 2009, MNRAS, 398, 2177 
Ma C. P., Maggiore M., Riotto A., Zhang J., 2010, MNRAS 

in press (arXiv: 1007.4201) 
Mo, H. J. & White, S.D.M., 1996, MNRAS, 282, 347 
Mori M., Burkert A., 2000, ApJ, 538, 559 
Park C, Choi Y.Y., Vogeley M.S., 2007, ApJ, 658, 898 
Sandvik H.B., Moller 

Sheth R.K., Mo H.J., Tormen G., 2001, MNRAS, 323 
Somerville R.S., Lemson G., Sigad Y., 2001, MNRAS, 320, 

289 
Springel V., White S.D.M. , Jenkins A., 2005, Nature, 435, 

629 
Oemler A., 1974, ApJ, 194, 1 
Padmanabhan N., Schlegel, D.J., Finkbeiner D. P. et al., 

2008, ApJ, 674, 1217 
Postman M., Geller M.J., 1984, ApJ, 281, 53 
Tinker J.L., Conroy C, Norberg P., 2008, ApJ, 686, 53 
York D.G., Adelman J., Anderson J.E., 2000, AJ, 120, 1579 



Wang H.Y., Mo H.J., Jing Y.P., 2007, MNRAS, 375, 633 
Wang L., Jing Y.P., 2010, MNRAS, 402, 1796 
Wechsler R.H., Zentner A.R., Bullock J.S. 
White S. D. M., 1996, in Schaeffer R., Silk J., Spiro, M., 

Zinn- Justin J., eds, Les Houches Session LX, Cosmology 

and Large Sclae Structure. Elsevier, Amsterdam, p. 77 
Weinmann S.M., van den Bosch E.G., Yang X.H., Mo H.J., 

2006, MNRAS, 366 2 
Zhang Y.C., Yang X.H., Faltenbacher A., 2009, ApJ, 706, 

747 

This paper has been typeset from a TgX/ lOT^X file prepared 
by the author. 



APPENDIX A: SCALE DEPENDENCE OF 
COLOR-DENSITY RELATION 

Previous studies have shown that age-related properties, 
such as galaxy color, are correlated with environmental den- 



sity and that this correlation is scale dependent ([Kauffmann 



et al.|2004| [Blanton fc Berlind |2007||Park et al.|2007| ). Here 
we examine the scale dependence of the color-density rela- 
tion quantitatively and in more detail. 

For each galaxy, we calculate its environmental den- 
sity at two different smoothing scales, denoted as 6r^ and 
^2- We then analyze the dependence of galaxy color on 
(fei, ^2)- The results are shown in Figure Al, where Ri = 
2h~^ Mpc is taken to be the reference scale. In each panel, 
the black dots show the values of [ln(l + Sr^), ln(l + Sr^)] 
for all the galaxies, and the red lines are [g — r) con- 
tours. The green line denotes our dividing line for high 
and low density regions, which is the —45"^ line passing 
through the point of (1.0, do), where do is the correspond- 
ing value of ln(l + ^^2) given ln(l + ^2Mpc/i-i) — l-O- I^ 
is seen that for R2 ^ 5/i~^ Mpc, the color-density depen- 
dence, {g - r) - {S2Mpch-^^^R2)^ is mainly on ^2Mpc^-i, 
and the additional dependence on 5r2 is very weak. For 
R2 — Sh~^ Mpc, a certain level of 5^2 -dependence is seen 
in low density regions. For R2 = lh~^ Mpc, on the other 
hand, the 5^2 -dependence appears in high density regions. 
These results clearly demonstrate the scale dependence of 
the color-density relation. There exists a special scale where 
the environmental density affects galaxy colors the most. 
The value of this 'optimal adaptive smoothing scale' is 
^ 2h~^ Mpc for Mr 20, and varies somewhat from high to 
low density regions. Our detailed analysis show that the 
variation is from ^ 1.5/i~^ Mpc in high density regions to 
^ 2.5h~^ Mpc in low density regions, consistent with that 
shown in Figure 8 for the scale dependence of the color- 
density correlation coefficient. 

Our further studies find that a better way to quan- 
tify the environmental effects on a galaxy is to use adaptive 
smoothing, for example distance to its nt/i-nearest neighbor, 
which takes into account the variations of the environmen- 
tal density naturally. We have tested the distributions for 
the nt/i-nearest neighbor distances, where n = 1,2,3,4 are 
considered. We notice that the distance to the Srd-nearest 
neighbor for Mr20 sample is peaked at 2h~^ Mpc with a 
range mainly between ^ lh~^ Mpc and 3h~^ Mpc, there- 
fore it represents the 'peak scale' suitably and adapt ively, 
and we propose it as the optimal adaptive smoothing scale 
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Figure Al. Scatter plots of [ln(l + (5i^^),ln(l + ^R^)]- In each 
panel, the red lines are galaxy color contours and the green line 
denotes the dividing line for high and low density regions. 



for the Mr 20 sample. The corresponding optimal adaptive 
smoothing scales for the Mr 19 and Mr 18 samples are re- 
spectively 6th- nearest and 8th- nearest neighbors. 
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